Similarity of percolation thresholds on the hep and fee 

lattices 

Christian D. Lorenz * Raechelle May, and Robert M. Ziff ^ 

Department of Chemical Engineering, University of Michigan 
Ann Arbor, MI 

February 1, 2008 



Abstract 

Extensive Monte-Carlo simulations were performed in order to determine the precise values of 
the critical thresholds for site = 0.199 255 5 ± 0.000 0010) and bond {p^!^^ = 0.120164 ± 

0.000 001 0) percolation on the hep lattice to compare with previous precise measuremens on the fee 
lattice. Also, exact enumeration of the hep and fee lattices was performed and yielded generating 
functions and series for the zeroth, first, and second moments of both lattices. When these series 
and the values of pc are compared to those for the fee lattice, it is apparent that the site percolation 
thresholds are different; however, the bond percolation thresholds are equal within error bars, and the 
series only differ slightly in the higher order terms, suggesting the actual values are very close to each 
other, if not identical. 

1 Introduction 

The percolation model is used to describe many problems that include a connectivity probabilty, particu- 
larly flow through porous media [[1], 0] . The three-dimensional lattices that are often used to model porous 
media, include simple cubic, body-centered cubic, face-centered cubic (fee) and hexagonal close-packed 
(hep). 
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Figure 1: Comparison of the layers that form the hep and fee struetures. Both of the struetures have the 
empty circles (A) as one layer. However, the hep lattice has a layer of the darkened circles (B) above and 
below that layer; whereas, the fee lattice has a layer of the darkened circles (B) above and a layer of the 
X's (C) below. 

In a recent paper, Tarasevich and van der Marck pointed out that the critical thresholds for site and 
bond percolation on the hep lattice were equal within the known error bars to the thresholds for the fee 
lattice. However, the values for the hep lattice (pjf = 0.1199 24±0.000 05 = 0.120 15±0.000 05 

[§) are not as precise as those found for the fee lattice (pg = 0.199 236 5 ± 0.000 001 [|], pf"^ = 
0.120 163 5 ± 0.000 001 [§]). This raises the very interesting question of whether or not the thresholds 
of these two lattices might in fact be identical, or at least the same to a very high precision. 

The similarity of the critical thresholds for the hep and fee lattices could be explained by the relative 
similarity of the hep and fee structures. In Figure [1], the hep and fee structures are shown in relation to one 
another. The two structures share the same first two layers (labeled A and B in Fig. |I]). The difference 
between the two structures occurs in the third layer. The third layer of the fee is a unique layer which 
fills the holes in the first layer which were not filled by the second layer (labeled C in Fig. [T]), but the 
third layer of the hep is the same as layer A [^. The structure of the hep crystal can then be summarized 
as A, B, A, B, A, B, . . . and the fee crystal is A, B, C, A, B,C, . . .. Therefore, we pondered whether it 
took higher precision to see this strucutral difference or if that difference had no impact on the percolation 
thresholds of the two lattices. 

In order to extend this study to higher precision we used a growth or epidemic analysis to determine 
precise values of the critical thresholds of the hep lattice, which were then compared to those of the 
fee lattice. We also carried out an exact enumeration study, because we couldn't find any such series 
expansions for the hep and fee lattices in literature; we could only find an exact enumeration study of a 
modified fee lattice [Q] . 
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Table 1: Unit vectors used to describe the neighbors in the fee and hep lattices. 



Lattice 


Vectors 




fee 


(1,1,0), (1,-1,0), (-1,-1,0),(-1, 
(-1,0,-1),(0, 1,1), (0,-1,1), (0,1, 


1,0),(1,0,1),(-1,0,1),(1,0,-1), 
-1), (0,-1,-1) 


hep (in the plane) 


(1,0,0), (1,1,0), (0,1,0), (-1,0,0), 


(-1,-1,0), (0,-1,0), 


(A to B) 


(0,0,1), (0,0,-1), (0,-1,1), (0,-1, 


-1), (1,0,1), (1,0,-1) 


(B to A) 


(0,0,1), (0,0,-1), (-1,0,1), (-1,0 


-1), (0,1,1), (0,1,-1) 



In the following sections, we report on the determination of the new values of for the hep lattice and 
the details of the exact enumeration. The results are summarized and discussed in the conclusion section. 

2 Percolation thresholds 

Precise values of the thresholds for bond and site percolation on the hep lattice were found using proce- 
dures similar to those outlined for site percolation in [|]] and for bond percolation in [^. A virtual lattice 
of 2048'^ sites was simulated, using the block-data method first described in [^. We distorted both lattices 
so that all sites fell on a simple cubic lattice. On these lattices, we grew individual clusters by a Leath-type 
algorithm which used the unit vectors shown in Table |T] for the hep and fee lattices. For the hep lattice, the 
unit vectors in the plane will always be the same, but the unit vectors needed to check the nearest neigh- 
bors above or below a certain site depends on which level the site is located {A or B). In Table |l], the unit 
vectors required to check the nearest neighbors in the plane and when going from layer Aio B and from 
to A are shown for the hep lattice. The critical thresholds were identified using an epidemic scaling 
analysis. In order to determine the critical thresholds at the reported precision, about 2 x 10^ clusters were 
generated utilizing about 10^'^ random numbers, which required a few weeks worth of computer time on 
ten workstations. 

The simulation was used to find the fraction of clusters P(s, p) that grew to a size greater than or equal 
to s sites. When p is near pc, one expects P(s, p) to behave as 

P(s, p) ~ As2-V((p - Vc)sn ^ As^-^[l + C{p - pjs'^ + . . .] (1) 

where r and a are universal exponents [jl], We assumed the values r = 2.189 and o = 0.445, which 
are consistent with other three-dimensional percolation studies [^, |^, [n|, |T^. Plots of s'^~'^P{s, p) versus 
s'^ for site and bond percolation on the hep lattice were used to find the value of the percolation threshold 
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Figure 2: Plot of s"" '^P{s,p) versus s*^ for bond percolation on the hep lattice. The curves plotted here 
represent p = 0.120 170 0,p = 0.120 163 5, and p = 0.120 160 (from top to bottom) respectively. 
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Figure 3: Plot of s"^ ^P(s,p) versus for site percolation on the hep lattice. The curves plotted here 
represent p = 0.199 260 0,p = 0.199 255 5, and p = 0.199 250 (from top to bottom) respectively. 
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which corresponds to horizontal behavior for large s. The results are plotted in Figures ^ and ^ and imply 
the following values for the critical thresholds for site (S) and bond (B) percolation: 

p^f = 0.199 255 5 ± 0.000 001 

p5 = 0-120 164 ±0.000 0010. (2) 



For the fee lattice, we previously found the values ^: 




0.199 236 5 ±0.000 0010 



0.120 163 5 ±0.000 0010. 



(3) 



The site thresholds for these two lattices differ by only 0.000 019, which is statistically significant being 
nearly 10 combined error bars apart. The bond thresholds, on the other hand, are identical within the error 
bars. 

3 Exact enumeration studies of the hep and fee lattices 

The similarity of the thresholds for the hep and fee lattices led us to also carry out an exact enumeration 
calculation, to see how the series of the two lattices compare. In exact enumeration, the problem is to find 
gst, the number of clusters containing s occupied sites and t vacant neighboring sites or bonds. Knowing 
gst, one can find the number of clusters (per site) containing s occupied sites by 



(4) 



the total number of clusters per site. 



N = M,{p) = Y. 



n3 = Y^ gstp'q* , 



(5) 



s 



s.t 



the percolation probability. 



P{q)=p-Mi{p)=p-Y^ 



sns=p-Y, SQstp'q^ , 



(6) 



s 



s,t 



(for small q), and the susceptibility. 



S{p) = M,{p) = Y: s^n, = s^Qstfq' , (7) 



(a) 



(b) 



Figure 4: A triangular cluster of three sites (gray colored circles) on the (a) hep and (b) fee lattice. In (a), 
the empty circles represent vacant sites in the plane of the cluster, the black circles represent vacant sites 
above and below the plane. Neighboring this cluster, there are 9 perimeter sites in the plane, 6 perimeter 
sites in the plane above and 6 in the plane below, yielding t = 9 + 6 + 6 = 21 total perimeter sites. In 
(b), the black circles represent the sites above the plane, and the X's represent the sites below the plane. 
This lattice has 9 perimeter sites in the plane, 6 in the plane above, and 7 in the plane below, yielding 
t = 9 + 6 + 7 = 22 total perimeter sites. 

where q = 1 — p. To calculate all these quantities, it is useful to construct the generating function 



To find gst, we developed a rather simple enumeration method based upon the cluster-growth algo- 
rithm, but using a deterministic sequence to decide whether each successive site (or bond) is occupied or 
vacant. After a cluster was finished, we stepped back to the last vacant site and made it occupied (if s 
was below the cutoff), or stepped back through the last group of occupied sites to first vacant one before 
it, and made that site occupied (if s was at the cutoff), and then returned to the growth algorithm again. 
In this way we went through a binary search of all possible growth scenarios. This method is similar to 
the algorithm described by Redner [[T3|]. We tested our algorithm with published results [|T^ in 2- and 
3-dimensions, and found agreement. Although slower than Merten's Fortran code [|T^], our program was 
easy to code and generalize for the different lattices and both site and bond percolation. We found the 
various series to s = 9 in a few days of computer time. 

For bond percolation, we consider that s represents the number of occupied sites in a cluster, irrespec- 
tive of the number of bonds that are needed to connect them, and t represents the number of vacant bonds. 
We thus define gst as the number of clusters containing s occupied sites and t vacant perimeter bonds. 



G{p,q) = Y.9stp"q' ■ 



(8) 



s,t 
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Note that no bond is placed in internal, redundant locations. A cluster containing s occupied sites and t 
vacant perimeter bonds has a weight gstP^^^(f, because only s — 1 occupied bonds are needed to connect 
the s sites. Therefore, in bond percolation, all definitions (Q)-([7|) should actually be divided by p on the 
right-hand-side, or equivalently those definitions actually give pMn{p). 

First consider the gst itself, which we represent by the generating function (||), where the coefficients 
of p^ are the perimeter polynomials. Up to order s = 5, the results are: 

= pq^^ + 6p2gi8 + p3(8^22 ^ ^2^23 + 30g24) + ^4^2^24 + 27q^<' + 48q^^ + 96g28 + 144g29 + 158g=^°) 
+p5(24g28 + 6g29 + I32g30 + 26^1 + A23q^^ +780q^^ + 119^4 + I212g35 + 846^^6) + . . . 

GJ^^ = ^ Qp2ql8 ^ p3 ^^21 ^ Q^22 ^ 13^23 ^ 33^24) ^ ^4 (§^25 ^ 21g26 + 43^27 ^ 93^28 ^ 153^29 ^ 153^30) 

+ 12g28 ^ 3gg29 ^ ^^4^30 ^ 333^31 ^ 357^32 ^ 33^^33 ^ 1143^34 ^ 1278g35 ^ 840g36) + . . . 

G^^^ = pgi2^gp2^22^p3(8^30^^g^3i^42g32) {2q^^ + 6q^'^ +36q^^ +84:q^^ +219q^'^ +249q^^ +326q^^) 
+p^{30q^^ + 120g^5 + 372^^6 ^ 792^47 ^ 1596^48 ^ 2328?^^ + 3576?^° + 3072g^i + 2739^^^) ^ _ _ ^ 

G^^P = pgi2^gp2^22^p3(8^30^^g^3i^42g32) +p4(2g36 + 6g37+ 36g38 +84^39 ^ 219g^°+249g^^ +326^^^^ 
+p5(g42^4g^3^3V^ + llV5 + 357g^6 + 780g^^ + 1611g^s + 2382g49^3513g5° + 3090g5i + 2739g^2)^___ 

For site percolation, the differences between the two lattices start to show up with s = 3. For example, 
the term for s = 3 and t = 21 occurs for the hep but not the fee lattice. This corresponds to the cluster of 
three sites in a triangle on the plane, as shown in Fig. ^. For bond percolation, the first difference does not 
occur until s = 5. 

First we compare the moments for site percolation, writing each moment for the two different lattices 
adjacent to each other for easy comparison: 

M({s = p - + 8p^ + - 6p^ + 30p^ - + 105p^ + 79p^ . . . 
Mq J = p - + 8p^ +p^- 6p^ + 30/ + 2p^ + 33j9^ + 513/ . . . 

M{^s/p = + 12g^^ - I2q^^ + 24g22 _ I2g23 + 53^24 _ ^gg^25 ^ 222g26 _ I40g27 ^ 252g28 + . . . 



Mf J/p = + i2gi8 _ I2gi9 + 3g2i + 12g22 + 6g23 + 30g24 _ ^39^25 ^ 73^26 ^ 4^^27 ^ 44^28 ^ _ _ _ 
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M/J^ = P + 12p2 + 8V + 50 V + 3012p5 + 17142^^ + 96228p^ + 532028p^ + 2918388p'' + . . . 



M^^^ = p + 12p2 + 8V + 50V + 30lV + 17148p^ + 96072p^ + 533286p^ + 2911166p^ + . . . 

For the first moment, we report Mi/p, which gives somewhat simpler expressions for the q series than 
Ml . Note that the series are identical up to order for Mq, q^^ (three terms) for Mi/p, and for M2. The 
coefficients for M2 differ a small amount between the two lattices for higher order. Note that the series for 
Mi{q) is actually given up to order 38 by the enumerations for s < 9. 

For bond percolation, the corresponding series are: 

pM/s = p - 6p^ + 8p^ + 33p^ + 132p^ + 554p^ + 2514p^ + 13152p'^ + . . . 
pM^Jf = p - 6p2 + 8p^ + 33p^ + 132p^ + 553p^ + 2526p^ + 131 16p^ + . . . 

M/'g = + i2g22 _ ^2g23 + 24g-'^° + 54g32 _ 204g33 ^ ^26g^^ + 8g'^^ + 96g-'^^ - ?,2q^^ + 276g^° - 
768g^i + 608g^2 _ 1800^43 + 3066g^^ - 1304^^^ + 360^^^ - 480^'^^ + lOSGg^^ - 3360^^^ + . . . 

Mfl = + i2g22 _ i2g23 + 24^3° + 54g32 _ 204g33 + i26g34 + 8g36 + 96g38 _ 32g39 ^ 276?^° - 
768g^i + 613g^2 _ igQOg^s + 3036g^^ - 1314^^^ + 450?^^ - 480g^^ + lOGlg^^ - 3480?^^ + . . . 

pMl'^B = p+12p2 + 132p3+1356p^+13344p^+127548p^ + 1194864p^+11033256p^+100692522p'^+. . . 

piW^Jf = p+12p2 + 132p3 + 1356p^+13344p^+127548p^ + 1194944p^+11033544p^+100697070p^+. . . 

Here we find an even closer agreement between the series of the two lattices than we found for site 
percolation. The series for both pMo and pM2 agree between the two lattices up to order 6 and then differ 
by a very small amount for the three orders beyond that, while the series for Mi agree up to order q^^ 
(first 12 terms). While it is impossible to prove or disprove that the thresholds for the fee and hep lattices 
are identical from these results, they clearly suggest that if the thresholds are indeed different, then they 
should be much closer for bond percolation than site percolation, as indeed we have found numerically. 
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4 Conclusions 



As a result of this work, we have shown that the critical thresholds for site percolation on the two lattices 
are definitely different. The value of Pc on the hep lattice (0. 199 255 5±0.000 001 0) is nearly ten combined 
error bars away from the previously reported value [|^] for the fee lattice (0.199 236 5 ± 0.000 0010). 
Also, the exact enumeration series for the two lattices begin to differ at relatively low order terms for site 
percolation. 

On the other hand, even at high precision, the critical thresholds for bond percolation on the hep and 
fee lattices have the same value, within the error bars (0.120 164 ± 0.000 001 0). Although the series for 
bond percolation on the two lattices do differ, the difference is incredibly small and does not occur until 
higher order terms. While this difference does not rigorously rule out equaltiy of the thresholds, we guess 
that the thresholds are in fact slightly different, but by an amount too small to be seen in our numerical 
simulations. 
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